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interferometers at the two LIGO sites, since they are sensitive in a different, 
higher, frequency band. Additionally, the variable orientation of the ALLEGRO 
detector provides a means to distinguish gravitational wave correlations from 
correlated environmental noise. We describe the analysis underway to set a 
limit on the strength of a stochastic background at frequencies near 900 Hz using 
ALLEGRO data and data from LIGO's E7 Engineering Run. 
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1. Introduction 



One of the gravitational wave (GW) sources targeted by the current generation of 
ground-based interferometric and resonant detectors is a stochastic background of 
gravitational waves (SBGW), produced by an unresolved superposition of signals of 
astrophysical or cosmological originm^JE]- Direct observational limits can be set on a 
SBGW by looking for correlations in the outputs of a pair of GW detectors. This has 
been done using two resonant-bar detectors U, two "prototype" interferometers 
and is currently being done with the two kilometer-scale LIGO interferometers [51 [7] 
(IFOs). In this paper we describe the first known effort to set a limit on a SBGW with 
correlations between an IFO and a bar. This pair of detectors-the LIGO Livingston 
Observatory (LL0)21 in Livingston, LA and the ALLEGRO resonant bar detector|U] 
in Baton Rouge, LA-is separated by only 40 km, the closest pair among the ten modern 
ground-based GW detector sites. JOI This makes it an attractive pair for probing the 
stochastic GW spectrum around 900 Hz. In addition, the ALLEGRO bar can be 
rotated, changing the response of the correlated data streams to stochastic GWs and 
thus providing a means to distinguish correlations due to a SBGW from those due to 
correlated environmental noise. 

In Sec. 12 we review the standard measure of stochastic background strength and 
the cross-correlation technique used to search for a SBGW. In Sec. 13 we describe 
the previous limits set on a SBGW with ground-based GW detectors. In Sec. ^ we 
describe the key features of the LLO- ALLEGRO correlation experiment in general. 
Section |S1 describes the data taken during LIGO's E7 engineering run and the analysis 
currently underway using those data. Finally, SecEldescribes the prospects for future 
correlation experiments using ALLEGRO and LIGO science data. 



2. Stochastic Background Measurements 

A SBGW is assumed for simplicity to be isotropic, unpolarized, Gaussian, and 
stationary. Subject to these assumptions, the stochastic GW background is completely 
described by its power spectrum. It is conventional to express this spectrum in terms 
of the GW contribution to the cosmological parameter = p/ Pait- 

r, f f\ 1 dpGW f dpQ^si 

Pcrit «ln/ Pcrit df 

Note that r^cwl/) has been constructed to be dimensionless, and represents the 
contribution to the overall ficw P^r logarithmic frequency interval. In particular, 
it is not equivalent to dflQ^^/df. Note also that since the critical density Pcrit, which 
is used in the normalization of f^Gw(/)j is proportional to the square of the Hubble 
constant Hg JT], it is convenient to work with /i^QQi^Gwl/); which is independent of 
the observationally determined value of /iioo = loo kni/s/ Mpc • 

The standard method to search for such a background is to cross-correlate the 
outputs of two gravitational wave detectors If the noise in the two detectors 
is uncorrelated, the only non-zero contribution to the average cross-correlation will 
come from the stochastic GW background. In the optimally-filtered cross-correlation 
method (described in more detail in |21 112L I13j). one calculates a cross-correlation 
statistic 

Y= f dtidt2hi{ti)Q{ti^t2)h2{t2)^ I dfh\{f)Q{f)h2{f) (2) 
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where hi^2{t) are the data streams from the two detectors, hi 2(f) are their Fourier 
transforms, and Q{ti — ^2) (with Fourier transform Q{f)) is a suitably-chosen optimal 
filter. It is sensitivej^l to backgrounds on the order of 

^"""^-{t I df-J^Sm^) ■ (3) 



PPlU)P2{f), 

This sensitivity improves with time and is limited by the power spectral densities 
Pi,2{f) of the noise in the two detectors. The factor 

712(/) = diab ^ / d^n e»2^/^^ A.7c pafc(s^) 

in the numerator of the integral is the overlap reduction function |14|. which describes 
the observing geometry. Here P°j'(f2) is a projector onto symmetric traceless tensors 
transverse to a direction Cl and df^ are the detector response tensors for the two 
detectors. These are the tensors with which the metric perturbation hab at the detector 
should be contracted to obtain the gravitational wave strain h = d'^^hab- If Uo and Va 
are unit vectors pointing in the directions of an IFO's two arms, its response tensor is 

dab='^{UaUb-VaVb) (5) 

while the response tensor for a resonant bar whose long axis is parallel to the unit 
vector Wa is 

dab = WaWb . (6) 

The overlap reduction function is equal to unity for the case of a pair of IFOs (or 
an IFO and a bar) at the same location with their arms aligned, and is suppressed 
as the detectors are rotated out of alignment or separated from one another. The 
frequency dependence comes about for the following reason: if the wavelength is 
comparable to or smaller than the separation between two detectors, the detectors 
will see different phases of the wave at the same time, and this phase difference 
will depend on the direction of propagation of the wave. Since the stochastic GW 
background is assumed to be isotropic, averaging over different propagation directions 
suppresses the sensitivity of a pair of detectors to high-frequency waves. For example, 
a wave whose wavelength is twice the distance between the two detectors will drive 
them 180° out of phase if it travels along the line separating them, but in phase if 
its direction of propagation is perpendicular to this line. Figure ^ shows the overlap 
reduction functions for several detector pairs of interest. 



3. Previous Results 



The current best upper limit on a SBGW from direct observation with GW detectors is 
/iioo^Gw(900Hz) < 60 0], set by correlating the resonant bar detectors Explorer [T^ 
(in Geneva, Switzerland) and Nautilus'161 (near Rome, Italy).:}: A broad-band limit 
of h^QQ^Gwif ) < 3 X 10^ was set using a pair of "prototype interferometers" 0, and 
more recently an analysis of E7 engineering data from the LLO and LHO sites [Hj set 
a hmit of /ifool7Gw(40~215 Hz) < 7.7 x 10''. Analysis 7 of LLO and LHO data from 
LIGO's first science run (SI) is expected to improve substantially upon this limit. 

I The same group had Dreviouslv|17| set a hmit of /i^gQnQw(900 Hz) < 100 using a single bar 
detector. 
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Figure 1. The overlap reduction function for LIGO Livingston Observatory 
(LLO) with ALLEGRO and with LIGO Hanford Observatory (LHO). The four 
LLO-ALLEGRO curves correspond to four orientations: "N18°W" is nearly 
parallel to the y-arm of LLO ( "anti-aligned" ); "N108°W" is nearly parallel 
to the x-arm of LLO ("aligned"); "N63°W" is halfway in between these two 
orientations (a "null alignment" midway between the two LLO arms); "N48°W" 
is an intermediate alignment in which ALLEGRO was also operated during the 
E7 run. 



More stringent upper limits can be set on astrophysical grounds. They are 
detailed elsewhere [21 El j but we mention the bound from big-bang nucleosynthe- 
sismE], which states that a cosmological SBGW is limited by 

r ^hl,,now{f) < 10-' . (7) 

J10-**Hz / 

This broad-band limit tells us that any cosmologically interesting SBGW almost 
certainly lies several orders of magnitude below the existing limits. 

4. ALLEGRO-LLO Correlations 

Due to the 3000 km distance between the two LIGO detectors, the overlap reduction 
function (see Fig. ^| limits the range of frequencies at which they are sensitive to 
a stochastic background. The separation distance is half of a GW wavelength at a 
frequency of 50 Hz, so the upper end of the frequency range is a few times this. (The 
lower end of the frequency range is set by the seismic noise in the two detectors.) 
The ALLEGRO bar detector is far closer to the LIGO Livingston site than the LIGO 
Hanford site is, with only about 40 km separating the two and a "half- wavelength" 
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frequency of 3750 Hz. Thus the observing geometry for ALLEGRO-LLO allows for 
observations of correlations out to much higher frequencies. On the other hand, the 
sensitivity of ALLEGRO (see Fig. O is concentrated in two narrow frequency bands 
in the vicinity of 900 Hz, so correlations between ALLEGRO and LIGO Livingston 
probe a different part of the frequency domain than correlations between the two 
LIGO detectors. The major challenge in detecting a SBGW is accounting for the 
effects of noise which may be correlated between the two detectors. Even in widely 
separated detectors such as LLO and LHO there may exist correlations (the most 
obvious example being the due to the power line frequency and its harmonics) which 
could mimic or mask a real SBGW. With detectors in much closer proximity, such 
as ALLEGRO and LLO, the problem could potentially be more severe. A method 
proposed by Lazzarini and Finn ^H] is to rotate the ALLEGRO bar about a vertical 
axis and measure the cross-correlation for different alignments of the ALLEGRO-LLO 
pair. Changing the alignment of the bar from one of the interferometer arms to the 
other changes the sign of the expected correlation due to a SBGW. Thus any correlated 
noise which is independent of alignment can be removed. Another possibility is to 
record data with the bar aligned at 45° from the arms where no correlation due to a 
SBGW is expected. This "off source" measurement gives an uncontaminated estimate 
of the cross-correlated noise. 

LLO- ALLEGRO Coincidence Operation 

The ALLEGRO detector took data during LIGO "E7" Engineering Run that took 
place from December 2001 to January 2002. In the following section we discuss 
the initial investigations underway. Following the the E7 run ALLEGRO went off- 
line for a major upgrade. Beginning in the summer of 2002 a new transducer jl9| 
was installed which should substantially improve the sensitive bandwidth of the bar 
detector. Figure [5] shows the projected improvement in bandwidth and sensitivity. 
The installation was ongoing during the time of the LIGO SI science run in August- 
September 2002. Following the installation, a series of cool-down attempts were made. 
A succession of repairs were done on the vacuum system and a successful cool-down 
was finally made starting in February 2003. This allowed data to be taken with 
ALLEGRO during approximately the second half of the LIGO S2 run. 

5. Investigations with LIGO Engineering Data 

From December 28, 2001 to January 14, 2002, LIGO held its seventh Engineering Run 
(E7), the last engineering run before the start of scientific operation in August 2002. 
The ALLEGRO detector was in operation for the end of that run. Table^summarizes 
the ALLEGRO data which overlap with science-quality LLO data. 

If no excess correlations are found, and we assume that all of the ALLEGRO data 
have a level of noise given by the solid curve in Fig. El and all of the LLO data are 
of a quality described by [21], © tells us to expect to set an upper limit of around 
1.2 X 10''.§ 

Analysis of the data described in TableQis currently in progress. In the remainder 
of this section, we outline some of the technical challenges involved in this analysis. 

§ This would actually be better than the limit of 7.7 X lO** set in which is understandable because 
it involves only one engineering-quality LIGO data stream, while the LLO-LHO analysis involves two. 
The sensitivity of at both LIGO sites improved markedly between E7 and the science runs. 1221 
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Figure 2. The figure shows the ALLEGRO strain noise curve for operation 
during E7 (solid) 20 as well as the projected improvement for operation with the 
new transducer during S2 (dashed) 1191 
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7(907 Hz) 


Clean LLO Segments 


N48°W 


-0.43 


134 X 15min =33:30 


N18°W 


-0.90 


82 X 15min =20:30 


N63°W 


-0.03 


110 X 15min =27:30 



Table 1. Summary of ALLEGRO data coincident with LLO operation during 
E7. The table lists ALLEGRO'S alignment in degrees West of North, the resulting 
value of the overlap reduction function 7(/) (see Fig.0 in ALLEGRO'S sensitivity 
band, and the number of 15-minute segments of clean LLO data which coincide 
with ALLEGRO data in each alignment. 



First, the nature of the data recorded is different: LIGO data consists of a real 
time series hL{t), sampled at 16384Hz, which is related by a linear transfer function 
to the gravitational wave strain in the detector. The Fourier transform of this 

data stream contains frequencies from —8912 Hz to 8912 Hz, although the negative 
frequency components are all determined by the positive frequency ones. Because the 
ALLEGRO detector's sensitivity is confined to a band with A/ < /, its data stream 
hA{t) is heterodyned (mixed with an oscillating signal at frequency fn = 907 Hz) 
and low-pass filtered (before it is digitized) and then sampled at the reduced rate 
of 250 Hz. The result is to produce a complex signal hf{t) = e'^'^'^~^"*+'>''^^hA{t) 
whose Fourier transform h^{f ) « h^ifH + /)e*^'^'^-* contains completely independent 
data at frequencies —125 Hz < / < 125 Hz, which are approximations of the 
Fourier transforms of the data in the original time series between (907 — 125) Hz 



Stochastic Background Search with ALLEGRO and LIGO 



7 



and (907 + 125) Hz. 

The initial approach we are taking to correlating these data streams with different 
sampling frequencies is to heterodyne and resample the LIGO data in software 
(producing an h^{t) — e^^'^^~^"*~^'^^^hL{t) whose Fourier transform is h^{f) « 
hhifn + /)e^2'"^^) and then correlate it with the ALLEGRO data as though they 
were any two heterodyned data streams. This requires a certain amount of care, 
because the phase 0l at the start of a stretch of data of the reference oscillator for the 
LIGO heterodyning is specified (because it's done in software), while the phase 4>a of 
the ALLEGRO reference oscillator is only known if it is experimentally determined. In 
principle, it is not difficult to keep track of this, but it is also possible to post-process 
the cross-correlation statistics to account for ignorance of the phase difference (j)A~(f>L 
between the two oscillators. [25] 

Additional technical challenges include converting ALLEGRO data into the frame 
format used by LIGO [211 so that it can be processed by the LIGO Data Analysis 
System[3S|. Fortunately, the LIGO stochastic search codes^l EHl have been 
designed with an eye to working on both heterodyned and non-heterodyned data, 
and the data structures describing site geometry || allow for either interferometric or 
resonant bar detector geometries. 

6. Future Outlook 

The upper limit we expect to set from E7 data is several orders of magnitude above the 
current direct observational limit@] at the frequencies in question. But the practical 
application of the LLO- ALLEGRO cross-correlation method to engineering data will 
pave the way for future such observations, whose sensitivity is expected to improve 
for several reasons: 

(i) Most substantially, the sensitivity of the LLO IFO has already improved 
markedlyj22| and will continue to improve as the instrument approaches its design 
sensitivity PT]. 

(ii) allegro's new transducer will increase the bandwidth of the bar detector (see 
the dashed curve in Fig. [SJ. 

(iii) Future observing runs will provide more coincident data than the 80 hours or so 
taken during E7. 

Assuming a full year of coincident data in optimal alignment, with LLO operating 
at design sensitivity and ALLEGRO operating at the sensitivity given by the dashed 
curve in Fig. [3 we estimate using (O that an upper limit could be set in the range 
890 </ < 930 of hlof^ncwif) < 10~^ 
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